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Abstract 

The over-damped motion of a Brownian particle in an asymmetric, bistable, fluctuating potential 
shows noise induced stability: For intermediate fluctuation rates the mean occupancy of minima with 
an energy above the absolute minimum is enhanced. The model works as a detector for potential fluc- 
tuations being not too fast and not too slow. This effect occurs due to the different time scales in the 
problem. We present a detailed analysis of this effect using the exact solution of the Fokker-Planck 
equation for a simple model. Further we show that for not too fast fluctuations the system can be well 
described by effective rate equations. The results of the rate equations agree quantitatively with the 
exact results. 

PACS-numbers: 05.40.-a, 02.50.Ey, 82.20.Fd, 87.16Ac 

Models of over-damped Brownian particles in potentials with one or more minima and barriers serve 
as paradigms for many relaxation processes in physical, chemical, and biological systems. The minima 
represent the stable or metastable states of the system. Transitions from one state to another are induced 
by the interaction with the environment. This interaction is typically described by a thermal white noise. 
The dynamics of the system is dominated by characteristic time scales which are given by the mean first 
passage times for the escape out of the minima of the potential. The most simple system in this class of 
models is the problem of diffusion over a single potential barrier, pioneered by Kramers [p. 

In many situations, the potential fluctuates due to some external fluctuations, chemical reactions, 
or oscillations. The most prominent model of this kind is a Brownian particle in a symmetric bistable 
potential, subject to a harmonic force. It serves as the standard model for stochastic resonance In 
many other applications one has to consider stochastic, correlated fluctuations of the potential. Doering 
and Gadoua [||] investigated the situation of a symmetric, bistable fluctuating potential. They found a 
local minimum in the mean first passage time as a function of the barrier fluctuation rate. This effect 
has been called resonant activation and has been studied in detail by various people [@]-[|9|]. In most of 
these papers either a symmetric bistable potential or the escape over a single barrier has been studied. 
Escape rates for general potentials and dichotomous as well as Gaussian fluctuations of the potential have 
been calculated by Pechukas and Hanggi ||]]. Their results support a simple, physical picture of activated 
processes with fluctuating barriers: If the potential fluctuates fast, the rate for transitions over the barrier is 
determined by the average barrier. If the potential fluctuations are slow (static limit), the slowest process 
determines the rate. In an intermediate regime the rate is given by the average rate, which is greater than 
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the rate for fast or slow fluctuations. This picture has already been suggested by Bier and Astumian [Q| 
on the basis of a simple model with a dichotomously fluctuating linear ramp. 

In many applications, one does not have a single barrier or a symmetric, bistable potential. In a more 
general situation the potential will have several minima of different depth. In equilibrium, the system rests 
most of the time in the absolute minimum of the potential. But due to potential fluctuations, the position 
of the absolute minimum may fluctuate. Typical, biological examples of such a situation are membrane 
proteins like a cell surface receptor or an ion channel. When a ligand binds to the receptor, it changes 
the potential energy of the receptor and induces a conformational change of the receptor molecule. If the 
transition from one to the other conformation and back is always more or less the same, a description of 
this transition by a single coordinate may be sufficient. Then it is possible to model the conformational 
changes by the motion of a particle in a fluctuating potential. 

The effect of a periodic electric field on membrane proteins has been investigated theoretically [10] 
and experimentally [11]. Astumian and Robertson [10] described such a system by a two state model 
with periodically modulated rates. The effect of a periodic modulation can be related to a stochastic, 
dichotomous modulation [12]. This clearly demonstrates the relevance of our results to such biologically 
motivated models. We will come back to this point at the end. 

The motion of the over-damped particle in a fluctuating potential can be described by a Langevin 
equation 

x = f(x,t)+Ut) (1) 

where / = — We are using units where the friction constant and kg are unity. / (and V) depend on t 
since the potential fluctuates. ^(/) is a thermal (white) noise, it satisfies = 0, (^(f)^(f')) = 2T8(t — t'). 
In this letter, we restrict ourselves to the discussion of potentials with two minima, separated by a barrier. 
The position of the minima is ±x OT and does not depend on t. The maximum of V(x) is located at x = 0. 
The fluctuation of the potential is mainly a fluctuation of the depth of the two minima. Let us consider 
first the most simple, non-trivial version of such a model. Let us assume that the potential fluctuation 
is a dichotomous process and that V(x,t) takes the two different values V+(x) and V_(x). Further we 
assume that the absolute minimum of V+(x) (V-(x)) is the right (left) minimum. Such a model contains 
various time scales: Four mean first passage times for the two minima of V+(x) and V_ (x), the intra- well 
relaxation times for V+(x) and V-(x), and the characteristic time scales for the fluctuation of the potential. 
The mean first passage times and the intra-well relaxation times are fixed by the form of the potential; the 
fluctuation of the potential is an external parameter that can be varied. In a biological model for a cell- 
surface receptor, as mentioned above, it is determined e.g. by the concentration of the signaling molecule. 
Let V{x) = ±(V + (x)+V_(x)), AV(x) = ±(V + (x)-V_(x)). Then V(x,t) = V(x) +z(t)AV(x) where z(t) 
is a random process that takes two values ±1. Its static distribution is qo(z) = P+§{z — 1) + P-${z + 1). 
Let x be the correlation time of this process, so that (z(t)z(t')) = (z) 2 + (1 - (z) 2 )exp(-f/x). Without 
loss of generality we restrict ourselves to p_ < 1/2. What does one expect for such a model? Let us first 
suppose that the temperature is such that the typical barrier heights of the system are a few T. If X is small 
compared to the intra-well relaxation times of the potential, the systems can be described by an effective 
static potential (V)(x) = V(x) + (z)AV(x). The stationary distribution is po(x) =Cexp(— (V) (x)/T). In 
the static limit, the stationary distribution is po(x) = p+exp(— V+(x)/T) + /?_exp(— V-(x)/T). Suppose 
that p + is close to unity. Then the average potential is approximately given by V+(x) and po(x) is ap- 
proximately the same for small or large X. In the following we will show that between these two extreme 
situations an interesting effect occurs: The mean occupancy of the minimum at — x m , which is the mini- 
mum that has most of the time the higher energy, may become very large. We will show that this effect is 
related to resonant activation. 
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Figure 1: The mean occupancy in the left minimum of a piecewise linear, dichotomously fluctuating 
potential as a function of X for various p + . The parameters are T = 0.3, p + = 0.5, 0.6, . . . 0.9 , f+{x) is 
piecewise constant and takes the values 10, 1, -1, 2, -2, -10; f-{x) takes the values 10, 2, -2, 0.5, -0.5, 
-10. The values of x where f±(x) jumps are -4, -2, 0, 2, 4. The dashed lines are the results from @. 

To calculate po (x) or dynamic quantities of the system one has to solve the Fokker-Planck equation 



for this model. Here we assumed that the potential fluctuations can be parameterized by a single stochastic 
variable z{t). p(x,z,t) is the joint probability density for the stochastic variables x and z, and M z is the 
generator of the stochastic process z(t). To obtain the stationary distribution po(x) = J dzp(x,z), it is 
sufficient to analyze the stationary Fokker-Planck equation. A standard way to solve this equation is 
to expand p(x,z) in the right eigenbasis of M z . If the potential is piecewise linear, one obtains a set of 
coupled differential equations with constant coefficients, which can be solved analytically. The remaining 
task is to satisfy the continuity conditions for p(x,z), which is a simple linear algebraic problem. For the 
case of a dichotomous process, the force is f(x,z) = f{x)+zAf{x) as discussed above. In Fig. 1 we show 
results for the probability n = J° 00 po(x)dx of the particle to sit in the left minimum of the fluctuating 
potential as a function of the correlation time x, and for various values of p + . The choice of the potential 
is arbitrary, similar results can be obtained for other potentials as well. The results for h show that the 
qualitative discussion for small and large x given above is valid. But for intermediate x, h is much larger 
than expected. The system is able to detect fluctuations that are not too fast or not too slow. Such 
fluctuations enhance the occupancy of the left minimum, although it is most of the time not the absolute 
minimum of the potential. We thus observe a noise induced stability for the state which has most of the 
time the higher energy, at least when p + > 1 /2. The results in Fig. 1 show that for large values of p + this 
effect is even stronger. 

What is the reason for this effect? How can it be described quantitatively and how does it depend on 
the potential? To answer these questions, let us go back to the general case @. If x is large compared to 
the intra- well relaxation times for the two minima of V (x,z), the dynamics of the system can be described 
by an effective rate equation for the probability of the particle to sit in the left minimum, n, or in the right 
minimum, 1 — n. The rate equation is given by 



dp(x,z,t 




f(x,z) -T— )p{x,z,t) +M z p(x,z,t) 



(2) 



dn 
dt 



n {z)n + r 2 {z)(l-n) 



(3) 
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where n,2(z) is the escape rate for the left or right minimum of V(x,z). If AVi^z) = V(0,z) — V(^fx m ,z) 
is the depth of the potential, r\p,{z) x exp(AVi^(z)/T). For (|J) we can again discuss the two limiting 
cases of large or small x. For large x, the average occupancy is n(i —K>a)=n co = J dzqo{z)n{z) where 
n(z) = r\{z)/{r\{z) +^2(2)). For small X, the particle feels average rates and the mean occupancy is 
«(x = 0) =7io = n"/(n+^) where 77 = / dzrj(z)qo(z)- For the results presented in Fig. 1 the potential has 
been chosen such that hq is larger than and also larger than the occupancy determined by the average 
potential (V) (x). This explains qualitatively the x-dependence of n in Fig. 1. Let us now compare the 
results of the rate equation (|3j> with the results of the Fokker-Planck equation quantitatively. To calculate 
the stationary probability n as a function of x from <j3j) , we use the Fokker-Planck equation for the density 
p(n,z,t), 

Mn,Z,t) = d + +M (4) 

at an 

The stationary probability ti = Jq dnnpo(n) can be obtained from the stationary distribution p(n,z)- It is 
possible to calculate the stationary distribution po(n) = J dzp(n,z) for a dichotomous process jl4|], 



p (n) = C(n -n)(n- n-) a -~ x (n+ - n)^ -1 (5) 

where 

(p+r++P-r-)(n ± -n Q ) 

«± — 7 \ , (o) 

Xr + r_ [n + — «_) 

n= r2+ ~ r2 -. (7) 
r + — r_ 

C is a normalization constant. We introduced r± = r\± + r2±. r,± are the two values for the fluctuating 
rates r,(z), and n± = r2±/r±. The rates are given by (see [g], where T and L are set to unity. L e ff occurs 
instead of L, since we do not have a linear ramp.) 

AV 2 L 2 

n± = -j±-*L( exp (AV i± /T)-AV i± /T - I)" 1 . (8) 

Po(n) vanishes outside the interval between «_ and n + as it should have been expected, n does not lie in 
this interval. For h one obtains 

n = no + («oo — no) (9) 

x + x 

where x = y- + 7^. This shows that for a dichotomous process one always has a monotonic behavior of n 
as a function of x and the characteristic time-scale for the transition from no to is given by x. For more 
general noise processes it is possible to calculate n as well. The calculation is much more involved, but 
the typical behavior of h is the same as for the dichotomous case []TJ]. In Figs. 1 and 2 we compare results 
of the effective rate equations with results of the fluctuating potential. The agreement is indeed excellent 
for sufficiently large x. The value of x where the transition occurs, and the value of the maximum of h 
agree well with the exact results. The agreement between the rate theory and the exact results becomes 
better for smaller temperatures (see Fig. 2), which is clear since the rate equations are only valid for 
low temperatures. The motion in the average potential, i.e. the behavior of the system for small x, 
cannot be described this way, since the rates in the average potential differ from the averaged rates. The 
validity of the two-state model breaks down when x becomes smaller than the intra-well relaxation times. 
Nevertheless, we are able to understand why the occupancy h in the minimum that has usually the higher 
energy has the features shown in Fig. 1. Using the average potential, the average rate, and the average 
occupancy we are able to calculate analytically the three values of h. The transition between these values 
occur at x-scales given by the intra-well relaxation time and by x. This also explains the results in Fig. 2, 
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Figure 2: n as a function of X for p + = 0.7 and T =0.1, 0.2, 0.3, 0.4, 0.5 (from top to bottom). The 
potential is the same as in Fig. 1. The dashed lines are the results from @. 



which shows how the effect depends on the temperature. For lower values of the temperature, x becomes 
larger, and the values of n change due to the dependence of the rates on T. 

Comparing the above results with calculations for the mean first passage time shows that the noise 
induced stability is related to resonant activation. To calculate the mean first passage time, one has to 
introduce an absorbing boundary at the maximum of the potential and has to solve the Fokker-Planck 
equation with this boundary. The mean first passage time depends on the initial condition p(x,z,0), but 
usually the relaxation within the potential well is fast compared to the mean first passage time and the 
dependence on the initial condition is weak. As initial condition we choose p(x,z,0) = 8(x — x,-)8(z — zt) 
with Xj = ±x m . Let x = be the absorbing boundary. The solution of the Fokker-Planck equation is 



denoted by p(x,z,t |x;,z,-,0). The mean first passage time is then given by Q13| ] 



T MFPT = y dz J dxpi(x,z\xi) 



where 



pi(*,z|*i 



t^-p(x,z,t\xi,Zi,0) 



(10) 



(11) 



The average is taken with respect to the stationary distribution of Zi- For a piecewise linear potential 
Pi(x,z|x,) can be calculated using the same methods as for p(x,z) described above. Bier and Astumian 
[||] calculated the mean first passage time for a linear ramp, which is similar to our situation. They showed 
that for not too small x the system is well described by simple rate equations, as in our case as well. Using 
the rate equations one obtains for the mean first passage times eq. (15)) 



T MFPT ( T ) = x 0' + ( x ~* ~~ T 0i) ; 



-fx 



(12) 



where Xq; = Tf x and Xoo/ = p+/r,-_ + p_/r,-+. The x-dependence of the mean first passage time is similar 
to the x-dependence of h in (0). The characteristic time x has the same form as Xoo;. 

As already pointed out, it is possible to extend our calculations to various noise processes. The main 
qualitative features of the system are the same. Our results show that an asymmetric, fluctuating, bistable 
system can detect fluctuations that are not too slow and not too fast. For such fluctuations the occupancy 
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of the state that has usually the higher energy is enhanced. The results show that this effect may be 
very large, depending on the parameters of the system. For the largest value of p + in Fig. 1, the mean 
occupancy in the left minimum is very small for slow and fast fluctuations, but reaches a large value for 
intermediate values of T. If one lowers the temperature or modifies the potential it is possible to obtain an 
even larger effect, as shown in Fig. 2. 

To some extend noise induced stability can be compared to noise enhanced stability first found nu- 
merically by Day an et al [15] and observed experimentally by Mantegna and Spagnolo [|l6|], but there are 



several differences. The effect called noise enhanced stability in [16] is observed in a periodically driven 
system with a single, metastable minimum. The system remains in the metastable minimum for some 
time given by the mean first passage time for the barrier, and the mean first passage time has a maximum 
at some noise intensity. This effect is related to stochastic resonance. In our case the potential fluctuates 
stochastically with some correlation time x and has two minima. The less stable minimum is the abso- 
lute minimum in some configurations of the potential, but most of the time this minimum is metastable. 
Nevertheless it can be highly occupied. 

As metioned above Astumian and Robertson investigated a two-state model with periodically mod- 
ulated rates to describe the effect of an oscillating electric field on membrane proteins. Their results are 
in qualitative agreement with our results for the model with dichotomously fluctuating rates. One should 
expect that our results for the motion of a particle in a fluctuating potential, described by a Fokker-Planck 
equation, are relevant for such biologically motivated models. This is important, because the descrip- 
tion by a Fokker-Planck equation is much more general. Furthermore, for large frequencies or small 
correlation times the systems feels an average potential that cannot be described by fluctuating rates. 
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